Annual Scientific Report 


/// f6A’ V,L 

/ /V - V Lp 

3 


on 

INFRARED SIGNATURE MASKING BY 
AIR PLASMA RADIATION 


Grant No. NAG 2-1079 


Prepared for 

NATIONAL AERONAUTICS AND SPACE 
ADMINISTRATION 


For the Period 

June 1, 1997 to July 31, 1998 


Submitted by 

C. H. Kruger, Principal Investigator 
C. O. Laux, Associate Investigator 


August 1998 

Mechanical Engineering Department 
Stanford University 
Stanford, California 94305 









NASA 1998 - Kruger/Laux 


Contents 

CONTENTS 2 

1. INTRODUCTION 3 

2. EXPERIMENT 5 

2. 1 . Experimental setup and conditions 5 

2.2. Spectral Calibration Procedure 6 

Deconvolution of room air absorption features from the calibration signal: 6 

1. RADIATION MODEL 10 

3. 1 . NO ROVIBRATIONAL BANDS (FUNDAMENTAL AND FIRST OVERTONE) 1 0 

3.2. C0 2 BANDS: V 3 ANTISYMMETRIC STRETCH AT 4.3 JiM AND(V|+V 3 ) BAND AT 2.7 |J.M 1 1 

3.3. OH ROVIBRATIONAL BANDS 1 1 

3.4. Water vapor absorption 1 1 

4. COMPARISON OF MEASURED AND COMPUTED SPECTRA 12 

4. 1 . Range 2.4^t.2 |im 1 2 

4.2. Range 4.1 -4.9 |xm 14 

4.3. Range 4.9- 5.6 (xm 15 

SUMMARY AND CONCLUSIONS 17 

5. PUBLICATIONS (JUNE 1997- JULY 1998) 18 

6. HONORS AND AWARDS 19 

7. PERSONNEL 19 

8. REFERENCES 20 


2 



NASA 1998 - Kruger/Laux 


1. Introduction 

This report describes progress during the second year of our research program on Infrared 
Signature Masking by Air Plasmas at Stanford University. This program is intended to 
investigate the masking of infrared signatures by the air plasma formed behind the bow 
shock of high velocity missiles. This research is currently directed by Professor Charles 
H. Kruger, with Dr. Christophe O. Laux as associate investigator. This program is 
supported by a grant from the Ballistic Missile Defense Organization (technical monitor: 
Dr. David Mann) and is currently funded through the National Aeronautics and Space 
Administration (NAG 2-1079) with Dr. Stephanie R. Langhoff from NAS A- Ames as 
technical monitor. One Ph.D. candidate is currently involved in this program. 

Our previous annual report described spectral measurements and modeling of the 
radiation emitted between 3.2 and 5.5 pm by an atmospheric pressure air plasma in 
chemical and thermal equilibrium at a temperature of approximately 3100 K. One of our 
goals was to examine the spectral emission of secondary species such as water vapor or 
carbon dioxide. The cold air stream injected in the plasma torch contained approximately 
330 parts per million of CO 2 , which is the natural C0 2 concentration in atmospheric air at 
room temperature, and a small amount of water vapor with an estimated mole fraction of 
3.8xl0' 4 . As can be seen from Figure 1, it was found that the measured spectrum 
exhibited intense spectral features due to the fundamental rovibrational bands of NO at 
4.9-5.5 pm and the V 3 band of C0 2 (antisymmetric stretch) at 4.2-4.8 pm. These 
observations confirmed the well-known fact that infrared signatures between 4.15 -5.5 
pm can be masked by radiative emission in the interceptor’s bow-shock. 

Figure 1 also suggested that the range 3.2 - 4.15 pm did not contain any significant 
emission features (lines or continuum) that could mask ER signatures. However, the 
signal-to-noise level, close to one in that range, precluded definite conclusions. Thus, in 
an effort to further investigate the spectral emission in the range of interest to signature 
masking problem, new measurements were made with a higher signal-to-noise ratio and 
an extended wavelength range. 
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Figure 1. Measured IR emission spectrum of -3100 K, atmospheric pressure air 
(see 1996-1997 annual report). 
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2. Experiment 

2.1. Experimental setup and conditions 

A full description of the experimental facility was given in the 1996-1997 Annual Report. 
For the measurements presented in this report, experimental conditions were essentially 
the same as in our previous work with the following difference. First, slightly wider 
monochromator slits were employed to improve the signal to noise ratio (the entrance and 
exit slits were set at 1.0 and 2.8 pm, respectively, vs. 0.75 and 2.8 pm previously). 
Second, a larger amount of H 2 O was injected in the torch (4.5xl0‘ 3 vs. 3.8X10" 4 H 2 O mole 
fraction) as OH fundamental bands are expected to radiate at wavelengths between 2.5 
and 4.15 pm. Third, the spectral range was extended to cover 2.4 - 5.5 pm (vs. 3. 2-5.5 
pm previously) so as to capture NO first overtone and OH fundamental rovibrational 
bands. The measured temperature profile, shown in Figure 2, was determined as in our 
earlier work from Abel-inverted profiles of absolute NO (1-0) bandhead intensities. The 
current temperature profile is approximately 300 K higher than in our earlier experiments 
(the maximum temperature is now -3400 K vs. -3100 K), a difference on the order of 
experimental uncertainties. For the present conditions, concentration profiles of major 
radiating species are shown in Figure 2. These concentrations were obtained by assuming 
chemical equilibrium at atmospheric pressure and at the measured temperatures since, as 
discussed in our previous annual report, the plasma produced by the torch under the 
present conditions was found to be close to Local Thermodynamic Equilibrium. 
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Figure 2. Temperature and concentration profiles (major radiating species) for the conditions of the 
spectrum presented in Figure 3 and subsequent figures of the present report. 
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2.2. Spectral Calibration Procedure 

As in our earlier experiments, the calibration was conducted using a radiance standard 
(Optronics Laboratories model OL550 tungsten strip) traceable to NIST standards. The 
calibration procedure is somewhat complicated by absorption features due to the water 
vapor and carbon dioxide present over the 6-meter optical path that separates the plasma 
source or calibration lamp from the detector. These two species absorb a significant 
fraction of emission between 1.5 and 2.9 pm and at wavelengths > 4.0 pm. Both the 
calibration and the plasma torch emission spectra are affected by room air absorption. It 
could be thought possible to recover the calibrated plasma torch emission spectrum by 
simply dividing the measured plasma torch signal by the calibration signal. Such a 
procedure would be flawed as absorption lines of room air may occur at different 
wavelengths than plasma emission lines. In other words, the absorption lines observed in 
the calibration spectrum may not absorb any plasma emission. The correct procedure is 
to deconvolve room air absorption from the calibration signal, and then to divide the 
plasma emission signal by the corrected (i.e. room air absorption corrected) calibration 
signal. With this procedure, the calibrated plasma spectrum represents the spectrum 
emitted by the plasma and absorbed by room air. Thus the absorption features of water 
vapor and carbon dioxide remain part of the experimental plasma emission spectrum. In 
order to model this spectrum, it is then necessary to determine the spectral emission of the 
plasma and then to solve the transport equation over a path of room air. The 
deconvolution procedure is presented next. 

Deconvolution of room air absorption features from the calibration signal: 

The raw calibration signal (PMT voltage) is shown in Figure 3. The absorption features 
of room air species (CO 2 and H 2 O) are indicated. Absorption by CO 2 at 4.3 pm can be 
removed by a straight line interpolation. This approximation is certainly valid because 
the OL550 tungsten filament emission and the grating response are smoothly varying in 
that spectral range. The other absorption features occur over wider wavelength ranges 
and cannot be deconvolved using such a simple procedure. For these features, our 
deconvolution procedure consists in dividing the calibration signal by the spectral 
transmission spectrum (convolved with the instrumental slit function) of room air. The 
transmission spectrum of room air was obtained with the HITRAN96 database 12 with the 
following parameters: room temperature of 298 K, path length of 6 meters, water vapor 
mole fraction of 0.013 and CO 2 mole fraction of 0.001 (1000 ppm). 
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Wavelength [urn] 

Figure 3. Raw calibration spectrum obtained with the OL550 calibration standard. 

Results of the deconvolution procedure around 1 .8 and 2.7 pm are presented in Figure 4. 
The red curve represents the corrected calibration using the absorption coefficient 
determined with HITRAN. As can be seen in Figure 4, the deconvolution procedure at 
1.8 pm produces a smoothly varying trace. At 2.7 pm, the recovered spectrum still shows 
residual H 2 O absorption lines which cannot be removed by varying the H 2 O mole 
fraction. We believe that these residual features could be due to inaccuracies in the 
HITRAN database or to uncertainties on the instrumental slit-function which was 
deduced from the entrance and exit slit-width readings and the theoretical reciprocal 
linear dispersion of the monochromator. In future work, the slit functions will be 
measured experimentally with a monochromatic light source to increase the accuracy of 
the deconvolution procedure. The recovered calibration spectrum between 2.5 pm and 
3.0 pm shows a cusp and an S-shape variation that is likely due to the spectral response of 
this ruled grating. Since the grating spectral response should be relatively smooth, we 
replaced the part of the red curve between 2.5 and 3.0 pm with a fourth order polynomial 
fit (blue curve). 
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Figure 4. Absorption-deconvolved calibration spectrum in range 1.7-3.0 pm (The absorption 
features are due to room air H 2 0 and C0 2 ). 


Results of the deconvolution procedure between 4.5 pm and 5.5 pm are presented in 
Figure 5. In this range, room air absorption is mostly due to H 2 O. As can be seen from 
the figure, the deconvoluted calibration signal exhibits very smooth variations in this 
spectral range. The small oscillations are due to an etaloning interference by the quartz 
window of the OL550 calibration lamp. These oscillations and residual variations due to 
inaccuracies on the recovery procedure were finally smoothed out using adjacent 
averaging (blue curve). 

Figure 6 shows the summary comparison between the raw calibration spectrum and the 
deconvolved calibration spectrum. 
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5. Absorption-deconvolved calibration spectrum in range 4.6-S.6 pjn (The absorption 
features are due to room air H^O). 



Figure 6. Calibration spectrum corrected for water vapor and carbon dioxide absorption. 
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1. Radiation Model 

The measured emission spectrum, corrected for the spectral response of the detection 
system and calibrated in absolute intensity using the procedure discussed in the foregoing 
section, is presented in Figure 7. As in our earlier set of experiments, both the 
fundamental bands of NO at ~5 pm, and the V 3 band of CO 2 at ~4.3 pm are observed. 
But in addition the region 2.5 — 4. 1 5 pm exhibits lines of OH fundamental and of the NO 
first overtone (Av=2). We summarize next our modeling efforts for these radiative 
transitions. 



Wavelength [pm] 

Figure 7. Measured IR emission spectrum 
(present work: P=1 atm, maximum plasma temperature = ~3400K). 

3.1. NO rovibrational bands (fundamental and first overtone) 

The fundamental (Av = 1) and first overtone (Av = 2) rovibrational bands of NO (X 2 n) 
are clearly seen in the experimental spectrum shown in Figure 7. A detailed, accurate 
model of fundamental and overtone bands of NO (X 2 n) was implemented in NEQAIR2- 
IR, as previously described in Ref. *. The code determines rotational line positions by 
diagonalizing the Hamiltonian published by Amiot . 2 For each vibrational band, 
vibrational dipole moments M v v defined as: 

(M vV ) 2 = (f 4V 0 r)D e (r)'Fv (r)drf , 

where D e stands for the electric dipole moment function (EDMF), were determined using 
the accurate ab initio EDMF of Langhoff et al7 The Honl-London factor expressions 
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(corresponding to Hund’s case a ) recommended by Spencer et al 4 were employed. The 
model provides very accurate line intensities and spectral positions, which is of particular 
importance for high-resolution spectroscopic diagnostics and for the accurate simulation 
of absorption of fundamental and overtone NO bands by atmospheric water vapor. 

3.2. CO 2 bands: Vj antisymmetric stretch at 4.3 pm and ( VJ+V 3 ) band at 2. 7 pm 

Computations of the CO 2 band spectrum at 2.7 and 4.3 pm were obtained with the 
correlated-k model 5 based on the parameters of EM2C Laboratory of the Ecole Centrale 
Paris. 6 ’ 7 The correlated-k model is a narrow-band model in which the actual spectrum is 
replaced on each narrow band by the reordered spectrum, so that the spectral integration 
is carried out using typically 10 points instead of several thousands. In the case of CO 2 
IR radiation, this model affords radiative intensity predictions within a few percent 
accuracy. 7 The parameters used for the simulations presented in this report are based on a 
16-point Gaussian quadrature and a spectral decomposition over intervals of width 25 
cm' 1 . 

3.3. OH rovibrational bands 

The fundamental (Av = 1) bands of OH (X 2 n) are clearly seen between 2.5 and 4 pm in 
the experimental spectrum shown in Figure 7. As for the NO IR bands, an accurate line- 
by-line model of fundamental (and overtone) bands of OH (X 2 n) was incorporated in 
NEQAIR2-IR. Rovibrational term energies and line positions for the 1-0 and 2-1 bands 
of this transition are determined by diagonalizing the Hamiltonian of Stark et al. 8 (with 
corrections of Levin et al. 9 ). In future work, we also intend to incorporate additional 
bands originating in higher vibrational levels, with line positions for these levels 
determined using the Hamiltonian of Coxon. I0 . Transition probabilities of the OH IR 
lines have a very strong dependence on centrifugal distortions of the vibrational 
potentials. As a result, the P and R branches show an anomalous distribution with intense 
P branches and very weak R branches. In the NEQAIR2-IR model, we utilized the P- and 
R-branch transition probabilities determined by Holtzclaw et al. 1 1 As for the case of NO 
IR bands, this model provides the highly accurate spectral positions required for reliable 
simulations of absorption by atmospheric water vapor. 

3.4. Water vapor absorption 

In this work, the HITRAN96 database 12 was used to determine the transmittance 
spectrum of room air, based on a water mole fraction of 0.014 and a carbon dioxide mole 
fraction of 0.001 over a 6-meter optical path (for more details, see Section 4.3). The 
procedure for taking into account absorption is as follows: first, the emission spectrum of 
the plasma is computed at high spectral resolution (10 points/angstrom, or approximately 
10 points per line). Then, attenuation of this spectrum as a result of water vapor and 
carbon dioxide absorption is determined with Beer’s law, using the line strengths and 
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spectral broadening coefficients of H2O and CO2 (HITRAN96 database). Finally, the 
resulting spectrum is convolved with the instrumental slit function. 

4. Comparison of measured and computed spectra 
4.1. Range 2. 4-4.2 fim 

The spectrum measured over the range 2.4— 4.2 pm is compared with the spectral 
predictions of NEQAIR2-IR in Figure 8. As can be seen, all spectral features are well 
reproduced by the simulations, except perhaps between 2.8 and 3.0 pm where the model 
slightly underpredicts the measurements (possibly because we have not yet incorporated 
the (3-2) band of OH in NEQAIR2-IR). 



Wavelength [pm] 


Figure 8. Measured emission spectrum in the range 2.4 - 4.4 pm and 
comparison with NEQAIR2-IR simulations. 


The contributions of the various radiating systems of importance in this spectral range are 
shown separately in Figure 9. It can be seen that the P-branch of OH extends into the 
CO2 V3 band at 4.3 pm. Thus, OH emission may pose a problem for the detection of IR 
signatures in this spectral range. 
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The mole fraction of OH was determined by matching the measured absolute intensities 
of rotational lines of the P-branch of OH. The OH concentration determined in this 
manner was then used to infer, using chemical equilibrium relations, the mole fraction of 
water in the air injected inside the plasma torch. As already mentioned in the 
introduction, the mole fraction of water injected in the torch was thus found to be 
approximately 4.5x10 This value is significantly lower than the mole fraction of water 
vapor in room air (0.013). This difference is not surprising as the air injected in the torch 
was prepared by compressing atmospheric air at an earlier time when the relative 
humidity may have been lower, and/or significant amount of water vapor may have been 
removed by the water trap mounted at the exit compressed air tank The presence of this 
water trap is unfortunately required for separate experiments that use the same 
compressed air tank. 



Wavelength [iim] 

Figure 9. Contributions of NO first overtone, OH fundamental, and CO* (V1+V3) bands to the total 
emission spectrum in the range 2.4 - 4.4 pm. Note the anomalously weak OH R branch. These 
simulations incorporate room air absorption over a 6-meter pathlength. 

Absorption by water vapor in the 6 meter optical path is particularly significant between 
2.5 and 2.95 pm, as can be seen in Figure 10 where spectral simulations obtained with 
and without water absorption are compared with the measured spectrum. This 
comparison underscores the importance of computing highly accurate positions for all 
emission and absorption spectral lines, in the present case those of the NO overtone, OH 
fundamental, and HiO absorption bands. 
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Wavelength [gm] 

Figure 10. Comparison between the measured emission spectrum and NEQAIR2-IR simulations with 
and without absorption by water vapor in the optical path. 


4.2. Range 4.1 - 4.9 nm 

Figure 1 1 compares the measured CO 2 V 3 band spectrum with the predictions of the 
correlated-k model -' 5 presented in Section 2.2. The model also accounts for absorption by 
room air CO 2 over the optical path separating the plasma from the detector. This 
absorption is clearly responsible for near extinction in the range 4. 2-4.3 pm of the 
emission from low-lying rotational levels of CO 2 . The lines appearing at both edges of 
the absorbed region correspond to “hot” CO 2 rotational lines. The model appears to 
overestimate the measured CO 2 band intensity by approximately 30%. This discrepancy 
may be due to the fact that the parameters used in the c-k model are only valid at 
temperatures below 2900 K, which is a temperature slightly lower than the maximum 
temperature of the plasma considered here (max. temp. = 3400 K). Another possible 
reason for the discrepancy may be that the air injected in the torch contained a lower 
concentration of C0 2 than the typical 330 ppm. Excellent agreement between the 
simulations and the experiment would be obtained by assuming a CO 2 concentration of 
~240 ppm. Yet, we believe that the intensity differences between the measured and 
predicted spectra are more likely due to our extension of the c-k model beyond its range 
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of validity. These relatively small differences in absolute intensities have little bearing on 
infrared signature masking issues. 
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Figure 1 1. CO 2 spectrum computed with the c-k model parameters of the EM2C Laboratory,**’^ and 
comparison with the measured spectrum. Note the effects of absorption by room air CO 2 in the 
range 4.2-4.4 pm. 


4.3. Range 4.9- 5.6 fjm 

Figure 12 compares measurements and modeling results for NO fundamental bands. Two 
simulated spectra are presented, with and without incorporating the effects of water vapor 
absorption over the 6 meter path of room air between the plasma and the detector. By 
matching the depth of the water absorption features, we were able to determine the mole 
fraction of water vapor in the room to be approximately 0.014 (this value is slightly 
different from the 0.013 mole fraction used for the calibration recovery, because 
experiment and calibration were not conducted the same day). The room temperature and 
humidity were monitored (Oregon Scientific Instruments) during the scans of NO 
fundamental bands. The measured values of 25°C and 42% relative humidity correspond 
to a mole fraction of H 2 O of 0.013 (Reynolds and Perkins, 13 p. 622), which is in excellent 
agreement with the mole fraction of 0.014 inferred from the absorption features of the 
measured spectrum. 


As can be seen in Figure 12, the predictions of the model (with H 2 O absorption) are in 
excellent agreement with the measured spectrum. It should be noted again that both the 
measurements and computations are on absolute intensity scales. 
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Figure 12. NO fundamental band spectrum computed with NEQAIR2-IR (with and without water 
vapor absorption), and comparison with experimental spectrum. Note that a (small) constant value 
of 0.8 |iW/(cm 2 .sr) was added to the simulated spectra shown in the figure in order to match the 
offset, possibly due to underlying continuum radiation, that is apparent at 4.92 |xm. 
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Summary and conclusions 


Detailed measurements and modeling of the spectral emission of an atmospheric pressure 
air plasma at temperatures up to ~3400 K have been made. The cold gas injected in the 
plasma torch contained an estimated mole fraction of water vapor of approximately 
4.5x10 and an estimated carbon dioxide mole fraction of approximately 3.3x10' . 
Under these conditions, the minimum level of air plasma emission is found to be between 
3.9 and 4.15 pm. Outside this narrow region, significant spectral emission is detected 
that can be attributed to the fundamental and overtone bands of NO and OH, and to the V 3 
and the (V1+V3) bands of CO 2 . Special attention was paid to the effects of ambiant air 
absorption in the optical path between the plasma and the detector. Excellent quantitative 
agreement is obtained between the measured and simulated spectra, which are both on 
absolute intensity scales, thus lending confidence in the radiation models incorporated 
into NEQAIR2-IR over the course of this research program. 



Wavelength [pm] 

Figure 13. Summary comparison between measurements and models of infrared emission by an air 
plasma at temperatures of -3400 K and atmospheric pressure. The simulations incorporate the effect 
of water vapor and carbon dioxide absorption over the 6-meter length of room air between the 
plasma and the detector 
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